


f 


j! S nN 
iby 4 ‘ ts Me , 









Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1975 


Hydrofoil simulation in six degrees of freedom. 


Ribeiro, Jose A. G. Salinas Moreira 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/20890 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
F (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist | | Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 


lil \ KNOX appointed -— and published -—- scholarly author. 


http://www.nps.edu/library 






LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 


HYDROFOIL SIMULATION IN SIX 
DEGREES OF FREEDOM 


Jose A. G. Salinas Moreira Ribeiro 





DUDLEY KNOX LIBRARY =e 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 93940 














NAVAL POSTGRADUATE SCHOOL 
honterey, California 





Pore oti SIMULATION iN 
Seo eGhkees UF FREEDOM 


by 


Jose A. G. Salinas Moreira Ribeiro 


December 1975 


| Thesis Advisor: G. J. Thaler | 


Approved for puolic reitease; distribution unlimited. 











T1/16063 





SECURITY CLASSIFICATION OF THIS PAGE (®hen Data Entered) 


REPORT DOCUMENTATION PAGE ee ed 


2. GOVT ACCESSION NO} 3. RECIPIENT'S CATALOG NUMGER 












. REPORT NUKGER 







5. TYPE OF REPORT & PERIOD COVERED 
Master's Thesis 


December 1975 


GC. PERFORMING ORG. REPORT NUMBER 


€. CONTRACT OR GRANT NUMBER(e) 


&. TITLE (and Subtitie) 


HYDROFOIL SIMULATION IN 
SIX DEGREES OF FREEDOM 













7. AUTHOR(e) 










Jose A. G Salinas Moreira Ribeiro 





AM ELEMENT, PAOJECT, TASK 
WORK UNIT NUMSERS 





9. PERFORMING ORGANIZATION NAME AND ADDRESS 


Naval Postgraduate School 






Monterey, California 93940 


yt, CONTROLLIPG OF FICE NAME AND ADOMESS 12. REPORT DATE 

Naval Postgraduate School Decenhe Q 
13. NUMBER OF PAGES 
Monterey, Califormia 93940 1g 
a. MONITORING AGENCY NAME & ADORESS(M! different trom Contrciling Ollice) | 18. SECURITY CLASS. (of thie report) 

Naval Postgraduate School Unclassified 
ff ° ° TION/ DOWNGRADING 
Monterey, California 93940 18a. UEC DOWNGRADING 


16. DISTRIBUTION STATEMENT (ol thie Report) 


Approved for public release. Distribution unlimited. 


17. DIST RIQUTION ST ATEMERT (of the ebetract entered in Block 20, it vifferent from Report) 


18. SUPPLEMENTARY NOTES 


19. KEY WORDS (Cortinue on reverce aide tf NéEC6GCETYT and identity ty biock number) 


Hydrofoil Freedom 











ABSTRACT (Continues on roverae side if nccoacery and tcentify by block number) 
The heeling behavior of a fully-submerged foil hydrofoil craft 
turning, has been studied by producing a digital commuter simulation 
progran, using a Digital Simulation Language on the digital comouter 
IBM system 360/67. The High Point PC(H)-1 is used as the model of this 


Simulation, 
LE a ET EE 











DD ony; 1473  EdITION OF 1 NOV 68 Is CMSOLETE 


(Page 1) f Song, SECURITY CLASSIFICATION OF THis PAGE (Wnor Data Entered) 





Hydrofoil Simulation in 


Six Degrees of Freedom 


by 


Jose A. G. Salinas Moreira Ribeiro 


‘7 
Lieutenant, Portuguese Navy 


Submitted in partial fulfillment of the 


requirements for the degree of 


De oeesneOt es oCLENCE IN ELECTRICAL ENGINEERENG 


from the 


NAVAL POSTGRADUATE SCHOOL 
December 1975 





DUDLEY KNOY § 


ABSTRACT 


The heeling behavior of a fully-submerged foil 
hydrofoil craft turning, has been studied by producing a 
digital ccmputer simulation program, uSing a Digital 
Simulation Language on the digital computer IBN systen 
360/67.The High Point PC(H)-1 is used as the model of this 


Samulation. 





TABLE OF CONTENTS 


ie. Piano DULL ht ==—===== _ aaa aes oa Sa oS So 12 
ihe MC TEC OIC FELC) OUP DEGI RT Sy SSS BS aaa coe ar i 
Peewee Eb Auer O Ly wn K ANGE NT -~—=—=————————-—-—-—-——— 13 

Pema Oe Oy ton Cool o> a == a ee a a ne ee = 15 

eeeeo inane. Ole Oly LON —— = — SS SR Sm ee Pave 

ee eee eon erey « DOCr EY “COMPONENTS —--—--~--------- 25 

Cee ictiteOn mata AcCk AND ANGLE OF SIDE SLIP----------- 28 
Pee Vio eONmOr eT OnCEhs AND MOMENT TERNS---—----—--— ~--30 
Peo como mole lh ON === ————— 2a 6 ee ee eS ee 32 
Bee SOIL UNTO ESE) CS OE EG Baa OS a at ms 32 

eine eer er ser oe One ©) (a ee a = = — — —— —— 36 
ere em) 0 On ee — Sis 

1. Vertical Gyros a 32 

Ce OOO ee = = = = = 35 

eG COV ONRO Ge eS a ee S15) 

yey [SUSE SS SUSU a la a 36 

IDs PeaeneCwetincURNiNG CuAnACTERESTICS STUDY------------- 38 
a Been ae ee ers 38 
ioe d— Un iets] ——————— ——— ee 33 
PRnrPeCea= lo eKNots,  nhudder to —Starboard----------- 40 

Se Dect=GOmKnOLcS, RuUddeL to POCt==-=---—=-------- 41 

fee eet OOS = = = 42 

5. Speed=50 Knots--~=--~-----~-~-~------------------- ~43 
Piette Comore Gon SLeP COMMAND ---=-3-------------- 70 

1. Speed=30 KnotS~--- 23-9 n re rn rn rel 
Pee t= sok) ots—— = ————{— == ———— ———-— = — TZ 

3. Speed=40 Knots<-<-<-<------------------ = ae 23 

Pe oa Oe oC a. = === — == — = — — — = — = = — Poa a tes 

Comer GULeomnien THE SPLITTING OF THE FORWARD FOIL---95 
arReemeurynt ton Or THE CORRECT ROLL ANGLEW------~----- oF 





V. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDIES-104 


APPENDIX A SUMMARY OF COMPLETE EQUATIONS----~-----~--- --105 
COMPUTER PROGRAMS--~--~<- ~~ 2-9-7 no nn enn 110 
LIST OF REFERENCES-----~-- a n= += -- -- = ---- -- + ------ ------ 126 
TNITIAL DISTRIBUTION LiST----------------~-------~--~---- 127 





Io 
CENT 
DiF 
pus 
LiF 


NET 


Feet 


Z. 
Feet 


Dimensionless 
Dimensionless 
Dimensionless 
Dimensionless 
Pounds 
Pounds 
Pounds 
Pounds 


Pounds 


Pounds 
Pounds 
Pounds 
Pounds 
Pounds 


Pounds 


TABLE OF SYMBOLS 


Area of foil 


Area of strut 


Center of gravity 


Drag coefficient of a foil 

Drag coefficient of a strut 

Drier cOcrtticirent 

Side force coefficient 

Craft's centrifugal force 

Drag force on a foil (water axes) 
Drag force on a Strut (water axes) 
Lift force on a foil(water axes) 
Perpendicular force to craft's 


deck 


Side force on a strut(water axes) 
Force in direction of body X-axis 
Force on a foii in body X-axis 
Force in direction of body Y-axis 
Force on a foil in body Y-axis 


Force on a strut in body Y-axis 





ZiF 


ZS 


XX 


1G) 4 


ZZ 


Jaa 


Ler 


Pounds 


Pounds 


Pounds 


2 
Feet/Second 

Z 
Feet/Second 


2 
Feet/Second 


Radians 
Feet 


2 
Slug>Foot 
2 
Suc hoot 


2 
Slug-Foot 


Foot- Pound 


Feet 


Feet 


Force in direction of body Z-axis 
Force on a foil in body Z-axis 


Force on a strut in body Z-axis 


Component of gravity in body 
X-axis 
Component of gravity in body 
fa cas 
component of gravity in body 


Z~axis 
Helm angle 
Height of Height Sensor above 


instantaneous water surface 
Moment of inertia about body 
Kod xa S 
Moment of inertia about body 
Y-axis 
Moment of inertia about body 


Z-axis 
Moment acting to produce roll 
about body X-axis 


Distance from CG. along body 


X-axis to a point of application 
of foil force 


Distance from CG. along body 


Vacs tO a polnt Of application 


of foil force 





a 





Zoe 


Z1iS 


Rgir 
oud T 


Taal 


ELS 


- 


Feet 


Feet 


Foot- Pound 


Slugs 


Foot-Pound 


Radians/Second 
Radians/Second 
Radians/Second 
Feet 

Radians 

Feet 


Feet 


Pounds 


Feet/Second 


Feet/Second 


Feet/Second 


Feet /Second 


Feet/Second 


Feet /Second 


Distance from CG. along body 


i=aers to a point of application 
of foil force 


Distance from CG. along body 


Z-axis to a point of application 
of strut force 

Moment acting to produce pitch 
about body Y-axis 

Mass of the craft 

Moment acting to produce yaw about 
body Z-axis 

Roll rate about body X-axis 
Pitch rate about body Y-axis 
Yaw rate about body Z-axis 
Crait's turning radius 

Rudder angle 


Submergence of a foil in earth 


axes 


Submergence in earth axes of the 


assumed point of application of 
strut force 


Thrust in direction of body X-axis 
Velocity in direction of body 
X-axis 


Velocity along earth X-axis 
Foil velocity in direction of body 


X-axis 
Velocity in direction of body 
Y~axis 


Velocity along earth Y-axis 


Foil velocity in direction of 


OS 





W Feet/Second 
W Feet/Second 
E 
We Feet /Second 
1 
xX Feet 
E 
sé Feet 
iD 
Zi Feet 
E 
a Radians 
oh 
B Radians 
1 
3 
p Slugs/Foot 
@ Radians 
Vv Radians 
6 Radians 
Subscripts 
Cc 
Pp 
S 
M 
F Ropers to £oi] 
S Refers to strut 
i 


body Y-axis 
Velocity in direction of body 
Z-axis 


Velocity along earth Z-axis 
Foil velocity in direction of body 


Z-axis 


Distance along earth X-axis from 


reference origin 


Distance along earth Y-axis from 


reference origin 


Distance along earth Z-axis fron 


reference origin 
Angle of attack of a foil 


Angle of side slip of a strut 


Water density 
Roll angle 
Yaw angle 


Pitch angle 


Refers to center foil or strut 


RefeoEsmeouport £O12t Or strut 


Refers to starboard foil or strut 


Refers to mid foil section 


When i appears in the subScript to a symbol, it 


indicates that the symbol is to be repeated with i 


SMecee omuenyemecplaced by €,2,5 aud M, as indicated. 


10 





ACKNOWLEDGEMENTS 


I wish to express my Sincere appreciation to Prof. 
George J. Thaler of the Naval Postgraduate School. His 
knowledge, dedication, patience, continuous advice and 
orientaticn provided the necessary help to the realization 
of this paper. 

To Lit. John J. Uhrin III, USN, for his help on the 
computer programming work. 

To my wife, Alice, who, while competing with the 
computer, due to my long periods of absence, gave me tne 


necessary moral support to finish this job. 


14 





Ie INTRODUCTION 


ee ae ae ee wes er eee ee ee 


This thesis deals with the yaw motion and turning 
Gharacteristics, in calm water, of the hydrofoil craft 
PC(H)-1 High Point. This will be the model used in the 
computer Simulation program for this study. 

In previous theses, the system was made conditionally 
stable in two modes, the pitch-heave mode and the roll mode. 
In this thesis, keeping the system stable and working in six 
degrees of freedom, its behavior will be analyzed in three 
modes, pitch mode, roll mode and yaw mode, when the 
Simulated model is turning by the effects of a hela step 
command, acting on the astern flaps. This analysis will 
also include the study of rudder effects, aft starboard and 
aft port flaps effects and the heeling behavior of the 
system. In a further investigation the forward flap will be 
Split in two parts and pitching and rolling effects will be 
analyzed. 

Since slight changes in the manual control of the 
system occur in such short time intervals and provoke 
tremendous effects, the system necessitates an automatic 
control that will provide the necessary outputs to operate 
the model according to the given inputs. This autcmatic 
control, already existent, will be improved and adapted to 
the new motion of the craft. 

Finally the favourable roll angle for a certain Speed 
and a certain helm step command will be calculated and 


compared with the actual roll angle of the craft. 
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Ae FOIL TYPZS AND FOIL ARRANGEMENT 


There are two principal foil types used in the 
memetrvccion Of hydrofoil crafts, the surface~piercing foil 


and the fully-submerged foil,as shown in Figure 1. 
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Surface-Piercing Foil Fully-Submerged foil 
FIGURE 1 


Ae Oh OROVOIL CRAFT 
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The surface-piercing foil provides a self-stabilization 
for the craft in pitch, height and roll. On the other hand 
the fully-submerged foil cannot adjust itself to a changing 
lift. This can only be achieved by a change in incidence 
angle of the foil flaps. The fully-submerged foil requires a 
more sophisticated senSing arrangement to control that 
incidence.But the fully-submerged rerevab Tb has a great 
advantage, when the craft iS operating in heavy séa.The 
surface-piercing foil responds to relatively small changes 
in the surface of the sea and gives a very rough ride.On the 
other hand, the fuliy-submerged foil is much less affected 
by the effects of the surface of the sea.This paper will 
only deal with the fully-submerged foil craft. 

The model for this simulation, the High Point 
PC(H) ~1,uses the foils in canard arrangement, which has a 
secondary foil as forward foil and a main foil as aft foil. 
Both foils are supported by vertical struts.Figure 2, 
important for the understanding of the digital computer 
Simulation progran, shows the configuration of the forward 
and aft foil.The drawing and the dimensions of the craft and 


the foils are shown in Figures 3A and 3B. 
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Pee COORDINATE SYSTEMS 


Three coordinate systems will be used in this study of 
the hydrofoil craft dynamics. They are : the Body Axes, the 
Earth Axes and the Water Axes. Each of these systems is a 
right hand orthogonal systen. 

The body axes coordinate system (X,Y,Z), shown in 
Figure 4, has ats origin at the center of gravity cf the 
craft and is fixed relative to the craft with the xX axis 
running longitudinally through the craft. This system is 
initially coincident with the origin of the earth axes 
system, but it moves with the craft as time progresses. X is 
positive forward, Y iS positive to starboard and Z is 


positive downward. 


ss POSITIVE p 


> Za (KLLiv To 78D} 
Z| X-AXIS 
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FIGURE 4 


ORIENTATION OF THE BODY AXES WITH RESPECT TO THE 
Cee OmD RE ONIONS OF THE POSITIVE VELOCITIES 
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The earth axes coordinate system (X , Y , 2) is fixed 
E E 


relative to the earth's surface. The origin may be at any 
convenient location, provided the body system origin is 


imatiaily at the same point. The xX axis lies in the 
; E 


horizontal plane ; it 1S initially coincident with, and 


positive in the same direction as, the body X axis. The Y 


lies in the horizontal plane and is positive to starboard 


when the observer is facing the positive X direction, The 
E 


Z axiS 1S normai to the horizontal plane and 1s positive 


J 


downward. The craft's motion relative to the calm water 
surface, (constraint of this paper), is described in terms 
of this coordinate systen. 

mumecenerai~ the body axes are displaced from the fixed 
earth axes by the three Euler angles, 9 (pitch angle), ® 
{roll angle) and wW (yaw angle). Figure 6 defines tke 
positive directions of the Euler angles. The transformations 
between earth and body axes are outlined in Ref.3, by 
transformation matrices, however, only the results will be 
considered here. 

The water axes coordinate system is aligned with the 
relative water velocity vector and resolves into lift, drag 
and side force directions. All hydrodynamic model test data 
is taken and presented as plots in terms of this coordinate 
system. 

The hydrodynamic force and moment data are obtained in 
a water axis system which is always oriented with respect to 
the relative water velocity. These data must be transformed 
into the body axes in order that their effect will he 
properly included in the equations of craft motions. 

In general, the water axis orientation is always 
changing with respect to the body axes. This leads to the 


meea for a transformation matrix, also outlined in Ref.3, to 


1) 





resolve the lift, drag and side force values in the _ body 
axes. The orthogonal water axis coordinate system is defined 
in Figure 5. The angles of rotation between water and body 


axes are a and fB8 and will be defined later. 
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ORIENTATION OF WATER AXES COORDINATE SYSTEM 


20 










KARTH X-Z PLANL 


FARTH X93" PLease 
(HORIZONTAL PLA Ne) 


POSITIVE DPCFOR 6=%=0) 


EARTH YoT PLANE 
(VERTICAL PLANE) 


BODY Z-AXIS 


VIEW LOOKING AFT. 


FIGURE & 


Dit Ove On POsiiIyE DIRECTIONS OF EULER ANGLE 


Ol 





Gee SOUATIONS OF NOTION 


The development of the six eguations of motion for a 
hydrofoil follows the same procedure as those for a 
displacement type huli. This development is clearly outlined 
in Ref.3 and therefore will not be done here. 


, The equations of motion are summarized in this form : 


Translation 


ae a See ae ee 


J = (1/m) PF - OW + RV (2-1) 
v = (1/m) FP + PH = RU (en 
= (1/n)F_ + QU - PV (2-3) 
Rotation 
P = 1/1 [LOR (1 Sey OR OP) 1 (2-4) 
XX BZ YY XZ, 
0 = 1/1 [ Dy il = < I o25 
pee Ct ag’ ad (2-9) 
R= 1/1 hs-p ~ <t ~(QR-P)T | DG 
/ ae Q ( ve ee (Q ) eo ( ) 


Investigations with the craft simulation show that the 
terms RV,PV, and PW are very small and may be neglected. 


Similarly, the tern [f is less than 10% of I and less 
XZ XX 
than 4% of I . Consequently, the terms containing I may 
ZZ XZ 


be neglected for this craft. 

Equations (2-1) through (2-6) are linear and rotational 
accelerations in body axes. 

Given the previous equations, the velocity terms in 
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body axes are obtained by integrating each equation one 


U = f U dt (2~7) 


V = fi dt (220) 
W = fi at (2-3) 
P - fi at (2-10) 
Q - fs at (aed) 
R - fi dt (2~12) 


im order to calculate the craft's position in earth 


time. 


axeS coordinates, these velocity outputs must be transformed 
from body axes coordinates into earth axes coordinates. The 


transformation of the linear equations become ;: 


U = Ucos@cos¥tV (cos¥sinésind -siny cos¢) 
E | 
+W¥ (cos¥sinécos¢tsinvsin ¢ ) (2-13) 
V = Ucosé@ésinvt+V (cos¥cos¢tsin ¥ sinésing) 
+W (sinvsin@cos¢-cos vy sin¢) (2-14) 
WH = -Usin@+Vcosésindt+tWcosécos ¢ (2-15) 
E 


Only the craft's velocity in the downward direction,2 - 
ate) 


is required to calculate its height above the water surface. 


That is, the vertical position in earth axes follows: 


Z i Ww Oat (2-15) 
E E 





The craft's position in earth axes coordinates is 


obtained from : 
/ U dat (2-17) 
iG E 


Vat oa he 
ts (2-18) 


The Euler angles, @, 98 and W are also of interest 


>< 
i 


< 
i 


and can be found by integration of Euler rates b, © see 
The rates in turn must be derived from the body angular 
rates. The Euler rates are not easily obtained because they 
occur about axes which are not orthogonal. This fact can be 
appreciated by referring to Figure 6 and recalling how the 
angles were defined. W is a rotation about the Z axis, 0 
is a rotation about the Y axis, and ® is a rotation about 


the X axis. They result in the following set of equations : 


¢ = p+Wsin 0 (2-19) 
6 = Qcos¢-Rsin > (2-20) 
¥ =(Qsind+Rcos¢) cos6+ (6-P) sing (2-21) 


The Euler angles are obtained by integration of these 
equations. 

Studies of this craft indicated that during abnormal 
Situations, roll angles of 15 degrees and larger may occur, 
but under normal operating conditions, roll angles will not 
exceed 10 degrees and pitch angles are not expected to reach 
3 degrees. Therefore, the following assumptions are 
justified ;: 


sin @ =0 (2-27) 


cos 6 =1 (2) 
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D. CALCULATION OF VELOCITY COMPONENTS 


The various foil velocity components must be determined 


at this point in order to calculate the angles of attack and 


Side slip. The total velocity components of each foil in 


body axes can be expressed in terms of the craft linear and 


angular velocities by considering 


represented in Figure 7 and Table 


Center Foil 


=e gh Se ge Pe ee a Ee ee —— ee oe «ee 


=e SS Se Se Se — 


i 


oa ew ag ee 
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XCF. 


the craft geometry as 


(2-24) 


(2-25) 


DADE 


(2-27) 


(228) 


(2-23) 


(22 35)) 


(2-31) 


(2-32) 


(2-33) 


(2-34) 





Week Tt 20 (2-35) 


@=Point of application of foil forces 


O=Point of application of strut forces 


L = L 
ACE AGS 
L = L = L = L = L 
XPF XPS eo XSS XMF 
L = |L = L = L 
wor YSS Neda eS 
L = L = L 
LPF ZSF ZMF 
L a end, are variables and in general will not be 
ZPS ZSS 
equal. 
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PIGURE 7 


DEFINITION OF SYMBOLOGY 


26 

















a: 


VSI I we de2e, |e 








ag 0 a a 0 Xr UOT IVOr paclyeO 
JUTOd eATIOSTIA 
ANZ, JNA, AWK, ANZ 5 INK y IWX 5 quswSsS Ttog DIK 
SSZ, SSA, SSX, SSZ SSA, SSX, -4nz49 pavogizis 
ISZ, ISK, ISK dSZ 5 ASK, ASX y {tog pavoqazis 
SdZ_ Sdi, SdX, SdZ 5 Sd&, Sdx , nig .A0g 
IdZ,, ddi,, ddX,, Adz , Adi y JdX y Ttoq 10g 
S0Z, SOA, SOX, SZ » SOK » SOX y qnizg ra4ue9 
JZ, FOX, JOX,5 IOZ y TOR y JOX 5 {TtToq 194Ua9 
STXY-Z STXV~-A STXV-X STXV~-Z STXY~A STXY~X 


SOWCu 





7O WOTsecopraey 
JO ZUTOg 


OY DSATIETIY uoTyvottddy ue O32 [ea T[Veavg a0 


WUTOG 03 OD WOOLY azvUutTps00g BUOLY a0l07 LOZ TOquAs 


I Slavs 





eT 





foe ONGLE OF ATTACK AND ANGLE OF SIDE SLIP 


The angle of attack of a foil is shown in Figure 8A and 
is defined as the angle whose tangent is the total relative 
velocity between the foil and water in the bedy Z direction, 
divided by the total relative velocity between the foil and 
water in the body X direction. 


ao = arc tan(W /U ) (2-36) 
a 


If the foil is not aligned with the body X axis, then 
the foil has a fixed incidence angle and the total angle of 


attack becomes : 


a =a + a (2-37) 
total fixed 


The side slip angle is shown in Figure 8B and can be 


wate 


2 
See ediG. tan [ /\O _+V o+h 
elt = 


represented by : 


Z 
. (2-38) 
a bog Fe rl 


During actual craft operations the angles will be less 
than 10 degrees, so that the inverse trignometric functions 
can be replaced by quotients. Except at low craft speeds, 
the forward velocity of the craft, U0, is much greater than 
the velocities, V and W. Applying these simplifications, the 


eguations become ; 


a (W uy 0) a. 
L 


a4 fF i-fixed 


B =v ys 


1 eee 
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Peep XAPANSLON OF FORCES AND MOMENT TERMS 


With the aid of Figure 7 and Table I,the forces and 
moments for the craft are easily written. It is important to 
note that when a numerical value is inserted for one of the 
alphabetical dimensions, the Sign of the value must be 
positive or negative according to positive and negative 
directions from the axis system origin at the C.G. One of 
the assumptions made upon commencing Sinulation was that the 
struts contributed with zero lift and that the foils 
contributed with zero side force. Also note that horizontal 
foils make no contribution to the body Y force and vertical 
struts do not contribute to the body Z2 force. Therefore, 
all body Y motion originates at the struts and all body Z 
motion originates at the foils. | 


The forces and moment equations become ; 


F =F +F +F +F +F +F +F +ng +7 (2-39) 
X XCF XPF XSF XNF xXCS XPS xSS X xX 

F =F +F tF + hg (2-40) 
faeerCS™ YPS YSS Y 

F =F +F +F + F +mg (2-41) 


mueacr ZPF ZSF ZUF Z 


Boe)  U(E GD TE LL) OUP aL) OCP LL) (2-42) 
Pere Zievesee 9Y 2 PS Y¥2SS ¥ 2 CS 


pee yp T(E EL) UCP aL) UF UL) UF OL) 
ere eceeer tee z Xk SE 4 X Hr xX Z CF 


+(F L) +(F L) +(F L) +(F L ) +(FL) 
Kaez oP MeceSEe <7 °NP. xX Z CS YA 18S) 
+(F L) #2 1 (2-43) 
PaeSS KOZ e 
N=(F L) +(F L ) +(FL) Za 
ex CS X PS eo SS 


All of these forces have a hydrodynamic origin except 


for the thrust and gravity terms. The thrust 1s always 
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associated with the body axis and requires no further 
expansion. The gravity terms can be Simplified using small 
angle approximations for the Euler angles. 

The standard expression for hydrodynamic forces in 
water axis coordinates is ; 


7 = guac i= Ss (2-45) 
aa ak 


Where L,D,S represent lift,drag and side force respectively. 
The lift, drag and side force coefficients vary as functions 
of angle of attack, angle of side slip, submergence, 
velocity, flap deflections, elevator deflections and rudder 
deflections. 

Ideally, mathematical expressions would have been 
developed to correctly depict the interrelation of all 
Variables which affect the force coefficients. However, 
except for the conditions when relationships are linear, 
hydrofoil technology has not advanced to a position which 
would yield such expressions. The difficulties which are 
encountered in deriving:-a mathematical expression include : 

1) The occurrence of cavitation and ventilation, which 
are not completely predictable. 

2) The even more unpredictable cessation of cavitation 
and ventilation. 

3) The nonlinearity oz the hydrodynamic coefficients. 

4) Lack of sufficient test data to completely describe 
all of the above. 
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Fee SLMUULATION OBJECTIVES 


The objective of this Simulation was to study the 
dynamic characteristics of motion along and about each bedv 
axis. That 18S, perturbations should be introduced ae) 
separately excite : 

1) Motion along the X-axis 

2) Motion about the X-axis 

3) Motion along the Y-axis 

4) Motion about the Y-axis 

5) Motion along the Z-axis 

6) Motion about the Z-axis 

The following assumptions were made prior to beginning 
Simulation : 

1) Craft equations of motion are valid only for the 

foilborne mode. 

2) Weight of the craft remains constant. 

3) Hydrodynamic coefficients are based on fully wetted 

surfaces, 1.e. no cavitation or ventilation. 

4) The craft may be operated in calm water or in a 


seaway. 
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B. COMPUTER SIMULATION 
The computer program used was based on the computer 
preograms of Ref.1 and Ref.2. 

As in the preceding References, before the main program 
could be assembled, Fortran Subroutines were used to find 
the hydrodynamic coefficients of foil lift and drag, and 
strut drag and side slip. The SETUP subroutine stored as 
data points the Simulation curves of Ref.4&. The subroutines 
INTERP and INTRP1 are interpolation subroutines to obtain 
the proper coefficient from the above simulation curves. For 
INTERP, values of angle of attack and forward flap angle 
were used to enter the curves and obtain a value of foil 
drag. Angle of side slip and submergence were used to enter 
INTR?E1 to obtain strut drag. Listings of the subroutines are 
Shown in the computer program section. The main program was 
then assembled using Figure 9 as a reference for data flow 
and is also listed in the computer program section. In one 
of the parts of this main program, the automatic control 
block, were introduced some nodifications, needed for the 
main objective of this paper, the study of the craft's 
turning characteristics. This automatic control blcsck is 


Still open for future studies. 
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Gam AUTOMATIC CONTROL 

Control signals originate at the pilot's controls, at 
the motion sensors and at the position transducer. These are 
combined and processed in the electronics to produce the 
signal for the servo-valve, The pilot will be able to move 
the helm, a lever, or a knob to introduce steering, altitude 
and attitude commands ; the motion sensors sense errors 
between the commanded and actual craft motions and _ produce 
electrical signals proportional to these errors ; the 
position transducer provides the feedback signal so that the 
control surface stops moving when it has reached the 
position corresponding to the processed error Signals. 

It 1S now necessary to consider how the control system 
components affect the equations of motion through the 
control surface deflections. In general, a signal originated 
at the pilot or at the motion sensors will be modified by 
the individual dynamics of the components in the signal path 
before the control surface motion actuaily occurs. Motion 
sensors to be discussed are vertical gyros, rate gyros, 
accelerometers and height Sensors. 

1. Vertical Gyros 

A vertical gyro here will produce one signal 
directly proportional to the pitch angle of the craft and 
one Signal directly proportional to the roll angle of the 
Grait. 

2. Rate Gyros 

A rate gyro will sense an angular rate about one of 
the three body axes. The dynamics of the gyro are such as to 
produce an output signal which is related to the actual 
angular rate by a second order differential equation. 

Outputs of accelerometers are also assumed to relate 
the craft motions by second order differential equations. 

To identify the accelerations sensed by the 
instruments, at will be assumed that there is one 


accelerometer at the C.G., mounted such as to fmeasure 








lateral acceleration, and one accelerometer above each foil 


mounted in the X-Y plane so aS to measure accelerations 
parallel to the Z-axis. 
The lateral acceleration at the C.G. is : 
a, Fm (3-1) 
the 
accelerations in the Z-direction at points above the center, 


Taking angular accelerations into account, 


port and starboard foil are : 


= (F Se = C Bee, 

oer ( Be. on Q ( ) 

a = (FP /m)t+L .Q-L .P (3~ 3) 
ZP Z, XPF YPF 

a= (F /m)t+L .OFL ct ee 

25 ( af XSF 2 yor? (30%) 


4. Height Senso 


22 


ta 


Ultrasonic and sonar types are either in use or 
expected to be used. For this study it will be assumed that 


the dynamic characteristics of the sensor, whatever the 


type, 
instantaneous height of the sensor above the water. 


will be such as to give a direct indication of the 
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ee ee ee eee — eee ee em re ee ee = ee oo ee 


fen UDURR EFFECTS 


Por the hydrofoil steering control, the struts can be 
viewed as vertical foils so tnat steering can he achieved by 
meede ing a Strut, by rotating a hinged trailing edge or by 
Memamingqg an extension of a strut below a foil, (spade or 
ventral rudder), aS in the case of the PC (H)-1. 

The rudder on the forward strut exerts an influence 
upon the forward strut side force coefficient. Based on 
Figure 10 the necessary changes were nade in the computer 
Simulation program. To document the rudder's intluence in 
the system and to show the craft's behavior, several figures 
are presented, corresponding to runs uSing different speeds 
eee successive rudder angies of 30,25,20,15,10, and 5 
deqrees at each speed. The craft started moving straight 
ahead, the rudder angie being introduced after 3 seconds, 
moo all the cases. 

1. Speed=30 knots 

For a forward speed U=30 knots, figure 11A shows the 
pitch angle of the craft, which had an oscillation at the 
beginning of the run, but it went to a steady state value of 
feo aedreesS, With the craft's bow up, in 4.5 seconds, for 
all the six different rudder angles used. The resuit was 
natural, because the rudder plays a negligible role in 
S@ramging the value of the foil lift coefficient and 
mreretore in the craft's pitch motion. 

For the same speed, in Figure 11B the different roll 
angies corresponding Bommiomibenredt |ruddersangies were 
represented. When the ship went straight ahead, naturally 


the roli angie was 0 degrees, but as soon as a rudder ale 
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of 30 degrees was introduced after 3 seconds, the craft 
started rolling, going to a steady state value of -1.29 
degrees for the roll angle in 19.5 seconds. On the other end 
of the scale for a rudder angle of only 5 degrees, the 
reaction was smaller, the initial oscillation was also 
smaller and the craft went to a steady state value of -0.216 
degrees for the roll angle in 13 seconds. Comparing with 
Figure 4, these negative values show a roll to the port side 
direction, although the ship was turning to starboard. This 
result would be perfectly acceptable for a conventional 
Ship, but banking in the opposite direction is desired for 
the hydrofoil because of the high speeds of operation. 

In Figure 11C were represented the several heading 
angles, with positive values, showing that the craft was in 
fact turning to starboard. At 3 seconds when 30 £4degrees 
rudder angle were introduced, the craft had a reaction, 
reaching a maximum heading angle of 2.02 degrees, but 
returning to a steady heading value of 1.39 degrees in 23 
seccnds. For a rudder angle of only 5 degrees the maximum 
heading of the craft was 0.33 degrees , getting a steady 
heading of 0.23 degrees in 17.5 seconds, because the rudder 
Change was smaller and then the reaction was smaller. 


. The craft's height above the water Z was 
E 


represented in Figure 11D. Starting with a value of 0 feet, 
the craft went to 3.9 feet above the water in 4.5 seconds. 
Comparing with the dimensions in Figure 3A, it was observed 
that the foils were well inside the water. Again it was 
observed that the rudder plavs a negligible role on the lift 
memron Of the craft, because the craft's height above the 
water was the same for all the rudder angles. 

The forward flap was represented in Figure 11E. it 
Started with a value of 0 degrees, reaching a peak of 13.9 
degrees to help the transient period. After 5 seconds it 
reached a steady state positive value of 10.8 degrees. This 


positive value means that the flap was oriented in the down 





direction and therefore offering more resistance to the 
water, provoking more lift on the bow part of the craft. It 
can be concluded that the main roles of the forward flap are 
craft's lifting and to keep a stable pitch angle. 

In Figure 11F the stern mid foil segment was 
represented. After a transient period with a maximum value 
of -1.15 degrees, it reached a vaiue of 2.62 degrees in 4.5 
seconds, aS in the preceding case it waS oriented in the 
down direction and again itS main role is to help the lift 
motion of the craft. 

1E5N Figure 11G the stern starboard flap was 
represented. Its value was 0 degrees, before the rudder was 
introduced after 3 sSeconds.Then for a rudder angle of 30 
degrees, the flap had a tranSient period with a maximun 
value of -2.33 degrees, going to a steady state value of 
~1.86 degrees in 16 seconds. For a rudder angle of only 5 
degrees,the maXimum negative value was only of -0.38 
degrees, going to a steady state value of -0.31 degrees in 
13 seconds. Having the starboard flap angles different 
values for different rudder angles, 1t was concluded that 
the stern starboard and port flaps have influence in the 
craft's turns, therefore a completely different conclusion 
in relation with the conclusion presented for the forward 
flap and mid aft flap. 

On the next groups of figures the same elements @as 
before will be represented. Similar conclusions can be 
applied, therefore only the data obtained will be reported. 

2. Speed=36 knots, rudder to starboard 

The pitch angle was represented in Figure 12A. It 
reached a steady state value of 0.01 degrees in 6.5 seconds. 

In Figure 12B were represented the roll angles for 
the different rudder angles, showing again the banking in 
the wrong direction. For a rudder angle of 30 degrees, the 
roll angle had a maximum negative value of -1.29 degrees, 
going to a steady state value of -1.12 degrees in 15 


seconds. For a rudder angle of 5 degrees, the maxinun 
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negative roll angle was of -0.21 degrees, going to a steady 
state value of -0.18 degrees in 12.5 seconds. 

The heading angles represented in Figure 12C reached 
higher values due to a faster speed, increasing gradually 
With tine, meaning that with the increase of speed the 
medier effect on the craft's turning was larger. After 40 
seconds of run time with 30 degrees of rudder, the heading 
angle reached 8.5 degrees and for 5 degrees of rudder, the 
heading angle reached 121.42 degrees, clearly larger values 
mean tOr the 30 Knots speed. 

In Figure 12D the height above the water was 
represented, with a steady state value of 1.65 feet in 3.5 
Seeonds, after starting on 0 feet. 

The forward flap angle was represented in Figure 
12E. On the transient period it reached a maximum value of 
9.21 degrees, going to a steady state value of 7.7 =degrees 
in 5.5 seconds. These values are smaller than for the spsed 
of 30 knots, because the thrust is bigger and it iS easier 
memethe craft to lift. 

In Figure 12F the stern mid foil segment was 
represented. It reached a maximum value of -1.37 degrees in 
the transient period, junping to a steady state value of 
1.74 degrees in 7.5 seconds. These values are also smaller 
than for the 30 knots speed py the same reasons of the 
preceding paragraph. 

In Figure 12G the different aft starboard flap 
angles were represented. For a rudder angle of 30 deqrees 
the flap angle reached a maximum value of -1.94 degrees, 
going to a steady state value of -1.6 degrees in 16 seconds. 
For a rudder angle of 5 degrees, when the rudder was 
introduced the maximum value was -0.324 degrees, reaching a 
eeeady State vaiue of -0.27 degrees in 13 seconds. 

3. Speed=36 knots, rudder to port 

With the same speed but with the rudder oriented to 

the opposite board, port, there were no changes in the pitch 


mode, the only changes were verified with the roll angle, 
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angle, heading angle and aft starboard flap angle, 
represented in Figures 13A, 13B and 13C. The correspending 
figures with the rudder to starboard were Figures 128, T2C 
and 12G and comparing these two groups of figures, it can be 
observed that the roll angle, heading angle and aft 
starboard flap angle have exactly the same values but with 
Opposite direction, because on the Y body axis direction the 
craft is perfectly symmetrical. 
4. Speed=40 knots 

In Figure 14A the several roll angles for different 
rudder angles were represented.When a rudder angle of 3C 
degrees was applied after 3 seconds the craft started 
banking in the wrong direction again, having a transient 
oscillation with a maximum value of -1.18 degrees. [It went 
to a steady state value of ~1.032 degrees in 15 seconds. 

In Figure 14B the 4@ifferent heading angles were 
represented. The craft went straight ahead and started 
turning when the rudder angle was introduced after 3 
seconds. The heading angles increaSed gradually with tine, 
for a rudder angle of 30 degrees the heading angle reached 
the value 11.9 degrees after 40 seconds of run time. For a 
rudder angle of 5 degrees the maximum heading angle reaciaed 
1.99 degrees in 40 seconds. 

In Figure 14C the height above the water was 
represented. It went to a steady state value of 0.628 feet 
above the water in 11 seconds. 

In Figure 14D the forward flap was represented. 
After a transient oscillation, with a maximum value of 6.9 
degrees, it dropped to a steady state value of 6.19 dgrees 
in 7 seconds. 

The aft starboard flap angles were represented in 
Figure 145. For 30 degrees rudder angle, the flap angle 
reached a maximum negative value of -1.78 degrees and went 
to a steady state value of -1.48 degrees in 16 seconds. for 
5 degrees rudder angle, the flap angle had a maximun 


negative value of -0.29 degrees and went to a steady state 
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value of -0.25 degrees in 11.5 seconds. 
5+ SEeed=290 knots 

1 Figure 15A the different roll angles were 
represented. For a 30 degrees rudder angle, the roll angle 
had a maximum value of -0.99 degrees and went to a steady 
State value of ~0.88 degrees in 10.5 seconds.For a 5 degrees 
rudder angle,the roll angle reached a maximum value of 
-0.165 degrees and went to a steady state value of -0.147 
degrees in 8 seconds. 

The heading angles were represented in Figure 15B. 
For a 30 degrees rudder, the maximum heading angle was 18.3 
degrees in 40 seconds of run time. For a 5 degrees tcudder, 
the maximum heading angle was 3 degrees after 40 seconds. 

In Figure 15C the forward flap angle was 
represented.It went to a steady state value of 3.76 degrees 
in 3.5 seconds. 

In Figure 15D the different aft flap angles were 
represented. For a 30 degrees rudder angle, the starboard 
flap angle reacned a maximum value of -1.51 degrees and went 
to a steady state value of -1.26 degrees in 14.5 seconds. 
For a 5 degrees rudder angle, the starboard flap angle 
reached a maximum value of -0.25 degrees and went to a 


steady state value of -0.19 degrees in 8 seconds. 
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Pee orfECTS OF A HELM STEP COMMAND 


Improving and completing the automatic control bleck of 
References 1 and 2 and based on the block diagram of Figure 
6.7 in Ref.5, a Simulated automatic control system was 
built, including the dynamic characteristics of the sensors, 
the dynamic characteristics of the compensation filters, 
inputs to the servo vaives and the dynamic characteristics 
of the hydraulic actuators. 

There wWwaS an important limitation and perhaps a certain 
inaccuracy in the values used as corner frequencies of the 
compensation filters and the values used aS gains for the 
control actuation servos, because References 1,2,3 and 4 
used the model PC(H)71 for their studies and Ref. 5 used a 
more recent model, the PHN. Of course there was some lack of 
information, but some interesting results and some data were 
obtained for future studies. 

As mentioned before,the main function of the forward 
flap was to lift the craft and keep it stable above the 
water in the pitch-heave mode. To confirm this, the heln 
doesn't actuate the forward flap, but only the starboard aft 
flap and the port aft flap, therefore these two last flaps 
will be responsible for the craft's turning in conjunction 
With the rudder and will also be responsible essentially for 
fmee craft's rolling and finally for the craft's pitching. 
The forward flap is only actuated according to its job by 
the forward accelerometer servo,by the height sensor servo 
and by the pitch angle servo. The two aft flaps are 
actuated by the respective accelerometer sensors, by the 
pitch angle sensor, by the roll angle sensor and by the 
manual helm control. The rudder servo is actuated by the 
roll angle servo, by a yaw rate gyro and by ae heading 
command disengaged or engaged. in this last case it will 
stop and turn back the rudder when the desired heading is 
reached. 


Several figures will be shown with the roll angles, 
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heading angles, rudder angles, starboard and port flaps 
angles for different speeds and under the actuation of helm 
step commands of 150,120,90,60 and 30 degrees to the 
starboard side, at each speed. 
1. Speed=30 knots 

In Figure 16A the roll angles corresponding to 
different helm angles were represented. The helm was noved 
to starboard after 12 seconds of straight ahead motion, then 
the craft started turning in the right direction, also to 
starboard. For a helm angie of 150 degrees the roll angle 
was 9 degrees after 24 seconds and for a helm angle of only 
30 degrees the roll angle was 3 degrees after 30 seconds. 

In Figure 16B the heading angles were represented. 
The craft started straight ahead and when a starboard helo 
WaS introduced, the craft started turning to starboard, 
reaching 21.2 degrees to starboard after 30 seconds for a 
helm angle of 150 degrees and only 3.4 degrees to starboard 
in the same amount of time for a helm angle of 30 degrees. 

Figure 16C shows the different rudder angles due to 
a helm step command. While the ship went straight ahead the 
rudder was also straight ahead, but when the helm was 
introduced after 12 seconds,the rudder turned to starboard 
on the craft's turning direction. Finally when the helm was 
removed after 28 seconds, the rudder went back to its normal 
position. For a helm angle of 150 degrees the rudder reached 
a value of 5.1 degrees to starboard and for a helm angle of 
30 degrees it only reached 0.78 degrees to starboard. 

Pienaguse oD and 162 the port aft flap and the 
starboard aft flap were respectively represented. When a 
helm angle was introduced after 12 seconds, the craft 
Started turning and banking. After those 12 seconds the port 
aft flap had positive values, (down direction) and the 
Starboard flap had negative values, (up direction). When the 
meteco ls turning tO Starboard, it must roll to starboard, 
therefore a lift on the port side must exist. The port flap 


being down, when the craft advances more force is exert2d on 





that flap by the water and there is a lift force on that 
Sidae. The starboard flap 1S up and on that Side there is a 
down force. When the helm goes to zero degrees after 28 
seconds, the craft rolls in the opposite direction and it 
can be verified by the figures that the flaps have also 
moved in oppcsite directions. 

2. speed=36 knots 

In Figure 17A the roll angles corresponding to 
different helm angles were represented. When the helm was 
introduced after 10 seconds, the craft started rolling to 
Starboard. For a helm angle of 150 degrees the roll angle 
reached 6.2 degrees in 18 seconds and for a helm angle of 39 
degrees the roll angle reached 2.55 degrees after 27 seconds 
and it started banking in the opposite direction when the 
helm went again to 0 degrees. 

The heading angles were represented in Figure 17B. 
When a helm angle of 150 degrees was introduced after 10 
seconds, the craft started turning to starboard, reaching a 
Maximum heading angle value of 37.3 degrees after 28 
seconds and for a heim angle of 30 degrees, the craft 
reached a maximum heading angle value of 6.9 degrees after 
28 seconds. 

The rudder angies were represented in Figure i17C. 
For a helm angle of 150 degrees, the rudder angle reached a 
value of 5.37 degrees to starboard after 24 seconds, then it 
Started going back in the opposite direction when the helm 
Was removed. For a helm angie of 30 degrees the maxinun 
rudder angle value was 1.44 degrees to starboard after 2! 
seconds, then, as before, it started returning in the 
Opposite direction. 

In Figure 17D and 17E the port and starboard flaps 
were represented. As before when the craft was banking to 
Starboard, the port flap had positive values (down 
direction) and the starboard flap negative values (up 
direction). Both flaps reversed the motions when the craft 


rolled back to its normal poSition after 24 seconds. 
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3. Speed=40 knots 

In Figure 18A the roll angles for different hela 
angles were represented. For a helm angle of 150 degrees, 
the roll angle reached a value of 9.38 degrees after 19 
seconds. For a helm angle of 30 degrees the roll angle 
reached a maximum value of 2.43 degrees after 23 seconds. 

In Figure 18B the different heading angles were 
represented. For a helm angle of 150 degrees the heading 
angle reached 26.3 degrees in 24& seconds and for a helm 
angle of 30 degrees the heading angle reached 12,12 degrees 
in 24 seconds. 

In Figure 18C the different rudder angles were 
represented. For a helm angle of 150 degrees the rudder 
angle had a maximum value of 5.83 degrees to starboard after 
21 seconds and for a helm angle of 30 degrees, the rudder 
angle had a maximum value of 1.16 degrees to starboard after 
21 seconds. 

In Figures 18D and 18E the port and starboard flaps 
were shown. Their movements were the same as before, 
according to the craft's rolling motion. 

4. sSpeed=50 knots 

In Figure 19A the roll angles for different hela 
angles were represented. For a helm angle of 150 degrees 
the roll angle had a value of 9.45 degrees after 15 seconds 
and for a helm angle of 30 degrees the roll angle had a 
maximum value of 2.45 degrees after 20 seconds, then it 
Started decreasing because the craft rolled in the opposite 
direction due to the removal of the helm command . 

In Figure 19B the heading angles were represented. 
For a helm angle of 150 degrees, tne heading angle reached a 
value of 38.7 degrees to starboard after 25 seconds. For a 
helm angle of 30 degrees the heading angle reached a value 
of 11.1 degrees to starboard after 25 seconds. 

In Figure 19C the rudder angies were represented. For 
a helm angle of 150 degrees, the rudder angle reached a 


maximum value of 7.78 degrees after 21 seconds and then it 
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started decreasing, moving in the opposite direction due to 
the removal of the helm command. For a heim angle of 30 
degrees, the rudder angle reached a value of 1.54 degrees 
after 21 seconds and again it started decreasing. 

In Figures 19D and 19E, the port and starboard flaps 
were represented. Their motions are the same as before, but 


With different values. 
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C. RESULTS WITH THE SPLITTING OF THE FORWARD FOIL 


The forward foil was split in two parts to analyze how 
it improves or prejudices the craft's pitching, the craft's 
Gelling and the craft's turning rate in the desired 
directions. 

Some changes in the computer simulaton program vere 
necessary to achieve this trial. Thus, instead of 
considering a unigue submergence for the forward foil, two 
different sukmergences, Secp and Secs, corresponding to the 
forward port foil and to the forward starboard foil were 
considered. These two foils were moved independently, thus, 
two different forward flap angles, Elvi and Elv3 were 
considered. Also two different forward Oa drag 
coefficients and two different forward fOlE ae 
coefficients were introduced. Finally the block diagram of 
Ref.5 was adapted to an extra forward foil, provoking a 
small change on the Simulated automatic control systen. 

With a helm step command to starboard, the craft 
Started turning to starboard and it aiso started rolling to 
starboard. As waS explained in part B of this chapter, the 
two forward flaps were moved in a manner Similar to that of 
mime art port flap and the aft starboard flap, thus, the 
forward port flap went down provoking an upward lift force 
on the port side and the forward starboard flap went up 
provoking a downward force on the starboard side. Several 
runs with different speeds and different helm angles were 
done. In a general way, the roll angle in the desired 
direction was improved by 10% and the turning rate was 
improved by 8%, but the pitch angle had values 3% sraller 
than the values obtained with a single forward foil, because 
[iene hnalt of the forward foil in the upward direction, there 
are two lift forces, (one positive and the other negative), 
prejyudicing the craft's bow lift. 

Finally, the forward flaps were moved in opfosite 


directions to those of the preceding paragraph, thus, the 


q 





hOrward port flap went up and the forward starboard flap 


went down.The results were worse than those of the preceding 


paragraph. 
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D. CALCULATION OF THE CORRECT ROLL ANGLE 


Many runs with the Clakt exposed to. several 
disturbances were made. Different roll angles according to 
the existent conditions were obtained, but now a question 
arises : are these roll angles sufficiently near the optinun 
value? 

A block called Calculation of the Correct Roll Angle 
was introduced in the computer Simulation program to answer 
this guestion. Knowing the craft's turning radius for a 
certain speed and a certain helm angle, the craft's 
centrifugal force can pe calculated. The craft's 
gravitational force in the Z direction was already known, 
then, the force perpendicular to the deck, called net force, 
carn be calculated, by the gravitational and centrifugal 
forces. The angle between the gravitational force and the 
net force perpendicular to the deck is the correct roll 
angle. 

In Figures 20A and 20B the craft's trajectories in the 
X-Y plane for a speed of 40 knots and helm step commands of 
20 and 60 degrees respectively were represented. The craft's 
turning radius was calculated at the point where the heading 
angle was reversed 180 degrees in relation with the initial 
heading angle at the beginning of the run, Since this will 
be approximately the steady state value. 

ieraguce 20C the craft*s turning path, for a helm 
command of 160 degrees, a speed of 50 knots (extreme 
conditions) and a run time of 156 secondS was represented. 
The turning path waS approximately a circle, with an 
approximate turning radius of 410 feet. 

The values of the actual roll angle matched closely 
With the values of the correct roll angle, for the steady 
state period, although there were differences during tne 
transient period. These values are represented in Table II. 

In Figure 20D the craft's turning path for a helm angle 


of 160 degrees and a speed of 30 knots was represented. fhe 
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turning radius was smaller, 


The values of the actual roll angle matched again 


the values 


period. These values are shown in Table III. 
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TABLE IIi1 


Roll angles for helm=160 degrees and speed=30 knots 


Tine Actual roll angle Correct roll angle 
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MormmcoNGC LUSLONGmanD RECOMMENDATIONS FOR FURTHER STUDIES 
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The system was made stable in six degrees of freedon, 
with the following limitation : the dimensions, the foil 
areas, the weight and the thrust of the model PC(H)-1 were 
used in xeferences 1 and 2 and in the computer Simulation 
program. Due to a lack of information, because only very 
general information wasS available, certain values of the 
automatic control system of the model PHM, instead of the 
PC(H)-1, had to be used. 

A study of the craft turning characteristics was nade 
in this thesis. The simulated automatic control system was 
improved and adapted to perform this study. 

Roll was controlled by differential displacements of 
the aft flaps, actuated by the helm control, while pitch was 
essentially controiled by the forward flap. 

The rudder has negative effects on the hydrofoil's 
rolling direction, which should be the same as the turning 
direction, due to the high speeds of operation. 

It was not worth while to split the forward foil. 

Recommendations for further studies include : 

1. Explore more deeply the advantages or the 
disadvantages One the splitting of the forward foil, 
combining it with the motions of the aft flaps and the 
rudder. 

Pees y tie response of the craft in a Sinusoidal sea 
and in a random sea. Design improved control systems for sea 
State operation. 

3. Simulate the transformation from hullborne operation 


womtole Orne operation. 
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APPENDIX A 


SUMMARY OF COMPLETE EQUATIONS © 


EQUATIONS FOR ACCELERATION IN BODY AXIS 


me 

U = mx RYN (Aaa) 

v = Sr +Pw-Ru USB 
i 4 

W = 4F,+QU-PV | (A~3) 
m Z 

Doe : | 

p= ee a7 yy? BFS 17] (A-4) 

* — Be ——_ — ra tt 


VELOCITY IN BODY AXES 


U = f Udt Cha) 
a Vat Can) 
W= JS Wdt (A-8) 
P= f Pdt | (A-9) 
meee Od GEO) 
R= f Rdt CAaasis) 


TRANSFORMATION OF VELOCITIES FROM BODY TO EARTH AXES 


U,=Ucosdcosy+V(cosysindsing-sinycos¢ )tW(cosysindcosdt+sinysing ) 
(A-12) 


Vp=Ucosésinyt+V(cosycosp+sinysindsing )+W(sinpsindcosd-cosypsing ) 
(A-13) 


W,,=-Usino+Vcosdsing+Wcosécos¢ CAqs2) 
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fork) VELOCITY COMPONENTS IN BODY AXES 


Center Foil: 


Up = UtLy op @ (A-33) 

eee es -zcr * ’xcr@ = 

Wo = W-Ly op @ (A-35) 
Mort Foil: 

Up = Utbopp* Q-Lypp’ R saglO) 

Vp = V-Lyppt Ptlypy' R (A-37) 

Wp = Why ppt P-Ly pp" Q (A-38) 
Eretmpoard Foil: 

Ug = Utby ge Ghy ag R UR8S) 

ome zsh} hyoy’ ® oe 

Wo = Wty opt P-Ly opt Q (A-41) 
Mid Foil: 

Uy = UtLy yn? Q (A-42) 

ee ssa ame KEES) 

Wye = WeLyapQ (A-44) 
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MRANSPORMATION OF WATER PARTICLE ORBITAL VELOCITY FROM EART! 


moO BODY AXES 


Uwy = Uw cosycosé+tv_..,sinvcosé-W.., sine (A-45) 

Vwy = Unwy (cosvsinésind-sinycos® )+V,., (sinysin§singt+cospcossé ) 
Hy cosOsing | (A-46 ) 

Way = Upwy (cosvsinecoso+sinpsing )+V,., (sinpsin@cos}—cossing ) 
+Wry Cost cos (A-47) 


Each of these equations must be repeated for each foil. 


ReeAlIVE VELOCITY COMPONENTS 


Gai Uy ~ Uw | | (A-48 ) 
Ver = Vy ey Vy (A--49 ) 
We, = W, - Wy, (A-50) 


One set of equations for each foil. 


ANGLES OF ATTACK AND SIDE SLIP 


x 





= i _ * 
ay arctan Wey/Upy 2 Oy FIXED (A-51) 
V 
Uns Ving har 


e 


Pmone set of equations for each foil. Angles of attack and 


Side slip are calculated in radians. 
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mOLAL RELATIVE VELOCITY AT A PARTICULAR FOIL OR STRUT 


V = yULtV TW 7 (A-53) 


meORODYNAMIC FORCES IN WATER AXES 
(Each foil equation adn strut equation must be repeated for 


foil and strut, respectively) 


foal Lift: 
= Any2 2 
Fy, SV ALCe (A-54) 
Forl Drag: 
FF... = #pV7A_C (A-55) 
DF if 19 : 
Bim Drag: 
ak ZA aon i 
Fos s0V ACh (A-56) 
Berut Side Force: 
a 4 Je — 


TRANSFORMATION OF FOIL FORCES FROM WATER AXES TO BODY AXES 


PYIF = “Fi ppcosacossth, jpsine (A-58 ) 

_ . : mee 
Fyrp = ~Fprpsin8 (A-59) 
YOIF = ~FhypSinacosB-F) ;pcosa (A-GO) 
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